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ABSTRACT 

We report the discovery of LEHPM 2-59 as the coolest extreme M subdwarf (esdM) found to date. 
Optical and near infrared spectroscopy demonstrate that this source is of later spectral type than the 
esdM7 APMPM 0559-2903, with the presence of strong alkali lines (including Rb I), VO absorption at 
7400 A and H2O absorption at 1.4 fim. Current optical classification schemes yield a spectral type of 
esdM8, making LEHPM 2-59 one of only two ultracool esdMs known. The substantial space velocity of 
this object {Vgaiactic ~ —180 km s^^) identifies it as a halo star. Spectral model fits to the optical and 
near infrared spectral data for this and four other late-type esdMs indicate that LEHPM 2-59 is the 
coolest esdM currently known, with Te/ / = 2800-3000 K and -1.5 < [M/H] < -2.0. Comparison of T^f / 
determinations for M dwarfs and esdMs based on spectral model fits from this study and the literature 
demonstrate a divergence in Te// scales beyond spectral types '-~^M5/esdM5, as large as 600-800 K by 
types M8/esdM8. While this divergence is likely an artifact of the underlying classification scheme, 
it may lead to systematic errors in the derived properties of intermediate metallicity subdwarfs. We 
comment on the future of ultracool subdwarf classification, and suggest several ideas for addressing 
shortcomings in current (largely extrapolated) schemes. 

Subject headings: stars: chemically peculiar — stars: individual (LEHPM 2-59, APMPM 0559-2903) 
— stars: low mass, brown dwarfs — subdwarfs 



1. INTRODUCTION 

Subdwarfs are metal deficient stars lying below the 
stellar main sequen ce in optical color-magnitude dia- 
grams (|KuiDeilll939tl . Low mass subdwarfs typically ex- 
hibit halo kinematics {{V) = -202 km s^^; Gizis 1997), 
and are presumably relics of the early Galaxy, with ages 
>10 Gyr. With their extremely long lifetimes (far in 
excess of the age of the Universe), low mass subdwarfs 
are important tracers of Galactic structure and chemical 
enrichment history, and are representatives of the first 
generations of star formation. 

The optical and near infrared spectra of the 
coolest known M-type subdwarfs, like their solar- 
metallicity dwarf counterparts, are dominated by molec- 
ular absorption, including bands of CO, TiO, AlH, 
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11976; Leggett et al. 2000). Variations in elemental abun- 
dances (i.e., metallicity) can modulate these molecular 
signatures appreciably, affecting both the total molec- 
ular opacity (and hence overall luminosity at a given 
mass) and relative band strengths through differential 
chemical abundance patterns and modified atmospheric 

^ Visiting Astronomer at the Infrared Telescope Facility, which is 
operated by the University of Hawaii under Cooperative Agreement 
NCC 5-538 with the National Aeronautics and Space Administra- 
tion, Office of Space Science, Planetary Astronomy Program. 



chemistry. At optical wavelengths, metallicity effects in 
subdwarf spectra are seen succinctly in the relative band- 
strengths of metal oxides and metal hydrides; the former 
are weaker and the latter stronger in low er meta llicity 
dwarfs (Bidclman & Smcthells 1976: iMould fc Hvlandl 
11976: ..CottrelI ,197&: ■Bessell.49821 Current optical clas- 
sification schemes for M subdwarfs are generally tied to 
the relative strengths of these bands. For example, the 
most widely used scheme, that defined by Gizis (199"^ 
hereafter G97), divides metal-poor M stars into subdwarf 
(sdM) and extreme subdwarf (esdM) classes based on the 
relative strengths of CaH and TiO bands in t h e 6300 -7200 
A spectral region. G97 and IGizis fc Beidi (17991 have 
determined mean metallicities of [Fe/H] = — 1.2±0.3 and 
[Fe/H] = — 2.0±0.5 for these two classes of metal-poor 
dwarfs. 

Halo stars exhibit large space velocities relative to 
the Sun; hence, they are efficiently detected in proper 
motion surveys. Indeed, most low mass subdwarfs 
now known were originally identified in blue photo- 
graphic plate proper motion surveys, particularly those 
by Luyten (e.g., LHS and NLTT catalogs; Luyten 
1979a, b; see also Bakos, Sahu & Nemeth 2002 and Salim 
& Gould 2003). With the digitization of the red op- 
tical R- and /-band p hotographic plate UK Schmidt 
SER,C and AAO (, Harlev & Dawd IT98H ICann onI (1981 
Morga,n et a,lJI1992r. ESO llWest fc SchusteHII 982: Wc^ 



198411 and Palomar (POSS-I, Abell et al. 1959; POSS- 



II, Reid et al. 1991) sky surveys, new proper mo- 
tion surveys have begun to identify even cooler ob- 
jects, which emit more of their light at wavelengths 
redward of the visual band. These surveys - includ- 
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ing the Autom ated Plate Measu ring machine Proper 
Motion survey l|Scholz et al J 1200(1 he reafter APMPM) , 
the Calan-ESO proper motion survey uiz et al."2001'), 
the SUPERBLINK su rvey (Lcoinc. Shara, & Rich 2002; 
iLgpine & Shara"'2005') and the SuperCOSMOS Sky Sur- 
vey ( Hamblv et al. 2001a b.;d. hereafter SSS) - and as- 
sociated follow-up programs^ have pushed s ubdwarf dis- 
coverie s to the end of the M spectral class l)Scholz et alJ 
l2004aD. Metal-poor an alogs to even cooler L dwarfs 
(jKirkpatrick et all Il999() have been identified in the 
SUPERBLINK survey (ILeoine. Rich, fc S haral I2003al 
hereaf t er LRS03) and s erendipitously ( Burg asser et al.l 
I2003al 'Burga ssei< 120041) in t he Two Micron All Sky 
Survey ( Skrut skie eraLll2006l hereafter 2MASS). These 
low temperature ultracool subdwarfs, encompassing 
objects with spectral types sdM7/esd M7 and later 
(jBurgasser. Kirkpatrick. fc Lepind f2005l see also Kirk- 
patrick, Henry & Irwin 1997) provide new challenges for 
atmospheric modeling and extend our sampling of the 
halo population down to and below the hydrogen burn- 
ing minimum mass (Burgasser et al. 2003a). 

While several ultracool subdwarfs are now known to 
exist, only one ultracool extreme su bdwarf has been 
found, the esdM7 APMPM 0559-2903 ijSchweitzer et all 
IT999|) . To seek out even cooler subdwarfs, we have ini- 
tiated a program to obtain spectra of red proper mo- 
tion stars selected fro m the Liverpool Edinburgh High 
Proper Mo tion survey (iPokornv. Jones fc Hamblvll2003l: 
iPokornv et al. 2004, hereafter LEHPM). One of these 
sources, LEHPM 2-59, appears to be a later-type esdM 
than APMPM 0559-2903 based on both near infrared 
and optical data, and we identify it as the latest-type 
and coolest esdM found to date. In §2 we summarize the 
selection of this source from the LEHPM survey. In §3 
we describe near infrared and optical spectroscopic ob- 
servations of this and four other late-type esdMs, and 
describe their overall characteristics. In §4 we analyze 
line strengths, radial velocities and spectral types for the 
optical data, and estimate the distance and kinematics of 
LEHPM 2-59. In §5 we use spectral model fits to the ob- 
served optical and near infrared data to derive TeffS and 
mctallicities for the latest-type esdMs. We discuss these 
results in §6, focusing on the temperature scale of esdMs 
and future revisions to existing classification schemes for 
ultracool subdwarfs. Results are summarized in §7. 

2. SELECTION OF LEHPM 2-59 

The LEHPM catalog is a subset of 11289 proper mo- 
tion stars from the SSS detected in i?-band ESO and 
UK Schmidt plates covering 7000 deg^ of the Southern 
sky {S < —20°). This area excludes regions close to the 
Galactic plane and those fields with epoch differences less 
than 3 yr (the mean epoch difference is 8.5 yr). Catalog 
sources were selected to have O'.'IS yr^^ < A* < 20'.'0 yr~^, 
fJ,/(Ji_i > 3 and 9 < i? < 19.5, and were cross-matched 
with photographic Bj and In plates (via the SSS) and 
the 2MASS point source catalog. Further details on the 
construction and completeness of the LEHPM catalog 
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are given in IPokornv eTal] ll200l . 

We selected ultracool dwarf and subdwarf candidates 
from the LEHPM catalog based on red optical/near- 
infrared color and J-band reduced proper motion (RPM; 
Hertzsprung 1905; Luyten 1922), Hj = J+5 log^g fi+5 = 
Mj + 5 log;^Q Vtan — 3.73. RPM provides a measure of an 
object's absolute brightness and tangential space velocity 
(Vtan) independent of its (generally unknown) distance, 
and therefore enables the identification of intrinsically 
faint and/or halo stars. Optical/near infrared RPM dia- 
grams have become increasingly popular in searc hes for 
ultracool nearby and halo stars (Sahm fc Gouldl [2002). 
as these low temperature sources emit more of their flux 
outside of the optical bands. Here, we focus on the J- 
band RPM since the spectral energy distributions of late- 
type dwarfs peak at these wavelengths. Figure ^ displays 
the Hj versus {Re so ~ J) diagram for the LEHPM cata- 
log. The main cluster of sources running diagonally down 
the middle of the diagram is composed primarily of main 
sequence dwarfs. The smaller cluster of sources offset 
below and to the left is composed of metal-poor (shifting 
their colors toward the blue), high velocity (increasing 
Hj) halo subdwarfs. The smallest grouping in the lower 
left corner of the diagram is primarily composed of low 
luminosity white dwarfs. 

To select for the latest-type sdMs and esdMs, we ap- 
plied the following selection criteria to the LEHPM cat- 
alog: 

• Detection in both Reso £ind J bands, and 

• {Reso - J) > 1-5 and Hj > 19.25 or 

• {Reso ~ J) > 3.5 and Hj > 24.5 - IMReso - J)- 

The color/RPM criteria are illustrated in Figure ^ Im- 
posing {Reso ^ J) ^ 1-5 eliminates contamination by 
the (relatively few) cool white dwarfs; the remaining cri- 
teria cordon off the late- type extensions of the dwarf and 
subdwarf tracks. A total of 50 sources were selected in 
this manner, including LEHPM 2-59. Of these, 14 have 
been previousl y observed by other programs, including 7 
M8-M9 dwarfs fGizis'2002':'Lodicu et al."2005'), the sdM7 
SSSPM J1930-4311 (Scholz et al. 2004a, a .k.a. LEHPM 
2-31) and the esdM6 SSSPM J0500-5406 llLodieu et alJ 
120051. a.k.a. LEHPM 1-3861). Results for our complete 
LEHPM sub-sample will be presented in a forthcoming 
publication. 

3. OBSERVATIONS 

3.1. Near Infrared Spectroscopy 

LEHPM 2- 59 and the three late- type esdMs LP 
589-7 (Gizis fc PeidI 119991 esdM5). LSR 0822+1700 
jLfehie. Shara. fc RichI 12004 es dM6.5) and APMPM 
0559-2903 ijSchweitzer et al.lll999l e sdM7) were each ob - 
served with the SpeX spectrograph ll Ravner et al.l l2003D. 
mounted on the 3m Infrared Telescope Facility, during 
three runs on 2004 March 10, 2004 September 5-9 and 
2005 December 31 (UT). A log of observations is provided 
in Table ^ Conditions during these runs ranged from 
poor (cloudy and high humidity) during 2004 March, to 
clear with light cirrus during the 2004 September and 
2005 December runs. Seeing was typically 0'.'5-0'.'9 on all 
nights. Spectral data were obtained using the prism dis- 
persed mode, which provides low resolution 0.7-2.5 fim 
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spectra in a single order. Using a O'.'S slit, we obtained 
data with spectral resolution A/AA « 150 and dispersion 
across the chip of 20-30 A pixel^^. The slit was oriented 
to the parallactic angle in all observations to reduce dif- 
ferential color refraction, and the telescope was guided 
by spillover light from the targets in the imaging chan- 
nel. Multiple exposures of 180 s each were obtained in an 
ABBA dither pattern along the slit. Nearby AO V stars 
were observed immediately after the target observations 
and at a similar airmass (Az < 0.1) for flux calibration 
and telluric corrections. Internal flat field and Ar arc 
lamps were also observed with each target for pixel re- 
sponse and wavelength calibration. 

All spectral data wer e reduced using the SpeXtoo l 
package, version 3.3 (Gus hing. Vacca. fc Ravnal l2004() 
using standard settings. First, spectral images were cor- 
rected for linearity, pair-wise subtracted, and divided 
by the corresponding median-combined flat field image. 
Spectra were then optimally extracted using default set- 
tings for aperture and background source regions, and 
wavelength calibration was determined from arc lamp 
and sky emission lines. Individual spectral observations 
for each science target and AO star were normalized and 
combined using a robust weighted mean with 8a out- 
lier rejection. Data for the AO stars were also corrected 
for broad-band shape variations before combining. Tel- 
luric and instrumental response corrections for the sci- 
ence data w ere determined using the method outlined in 
iVacca et al . (2003), with line shape kernels derived from 
arc lines. Adjustments were made to the telluric spectra 
to compensate for differing H I line strengths in the ob- 
served AO V spectrum and pseudo- velocity shifts. Final 
calibration was made by multiplying the observed target 
spectrum by the telluric correction spectrum, which in- 
cludes instrumental response correction through the ratio 
of the observed AO V spectrum to a scaled, shifted and 
deconvolved Kurucz"^ model spectrum of Vega. 

The reduced near infrared spectra of the four esdMs are 
shown in Figure[21 These spectra exhibit the characteris- 
tic signatures of late-type esdMs, including blue near in- 
frared spectral slopes due to strong collision-induced II2 
absorption; shallow II2O absorption at 1.4 and 1.9 ^m; 
FeH and CrH absorption at 0.86, 0.99 and 1.2 /itm; weak 
bands of TiO at 0. 84 ^m; and K I and Na I fines at 0.77, 
0.82 and 1.17 /zm ijMould fc Hvland"1976t iLeggett et all 
I2OOOI) . Notably absent are the deep II2O and 2.3 /xm 
CO bands that characterize solar-metallic i ty late-type 
M dw arfs l|Baldwin. Frogel fc Perssoiilll973t Ijones et al.l 
11994 cf. Fig. 1 of Burgasser et al. 2004). The weakness of 
these bands makes near infrared classification of esdMs 
difficult (Burgasser et al., in prep.; see §6.2); however, 
the overall similarity of the spectrum of LEHPM 2-59 to 
those of LSR 0822-^1700 and APMPM 0559-2903, cou- 
pled with its stronger H2O absorption and bluer spectral 
slope, strongly suggests a late esdM spectral type. 

3.2. Optical Spectroscopy 

We obtained optical spectroscopy for LEHPM 2-59 and 
the three late-type esdMs LP 589-7, SSSPM 0500-5406 
and APMPM 0559-2903 on 2005 December 4 (UT) us- 
ing the Low Dispersion Survey Spectrograph (LDSS-3) 
mounted on the Magellan 6.5m Clay Telescope. A log of 
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observations is given in Table [21 LDSS-3 is an imaging 
spectrog raph, upgraded by M. G ladders from the original 
LDSS-2 (lAUington-S mith et alJlT99^ for improved red 
sensitivity. The instrument is composed of an STA0500A 
4Kx4K CCD camera that re- images an 8.'3 diameter field 
of view at a pixel scale of 0'.'189. A set of slit masks and 
grisms allow spectral observations at various resolutions 
across the optical and red optical band. For our observa- 
tions, we employed the VPH-red grism (660 lines/mm) 
with a 0'.'75 (4 pixels) wide longslit to obtain 6050-10500 
A spectra across the entire chip with an average reso- 
lution A/AA « 1800. Dispersion along the chip was 1.2 
A/pixel. TheOG5901on gpass filter was used to eliminate 
second order light shortward of 6000 A. Two long expo- 
sures were obtained for each target, followed immediately 
by a series of NeArHe arc lamp and fiat-field quartz lamp 
exposures. We also observed a nearby G2-G3 V star after 
each esdM target for telluric absorption correction. 

Data were reduced in the IRAF'' environment. We 
focus our analysis on the blue side of LDSS-3 (in two- 
amplifier readout mode), covering the 6050-8400 A re- 
gion. Raw science images were trimmed and subtracted 
by a median combined set of bias frames taken during 
the afternoon. The resulting images were divided by the 
corresponding normalized, median-combined and bias- 
subtracted set of flat fleld frames. Spectra were then 
extracted using the APALL task with background sub- 
traction but without variance weighting (i.e., not "opti- 
mal extraction"). The dispersion solution for each target 
was determined using the tasks REFSPEC, IDENTIFY 
and DISPCOR, and arc lamp spectra extracted using 
the same dispersion trace; solutions were typically accu- 
rate to 0.05 pixels, or 0.07 A. Flux calibration was de- 
termined using the tasks STANDARD and SENSFUNC 
with observations o f the spectral standard Hiltner 600 
l|Hamuv et alJll994 a.k.a. HD 289002) obtained on 2005 
December 3 (UT) with the same slit and grism combi- 
nation as the science data. Corrections to telluric O2 
(6850-6900 A B-band, 7575-7700 A A-band) and H2O 
(7150-7300, 8150-8350 A) absorption for each esdM/G 
star pair were determined by linearly interpolating over 
these features in the G star spectrum, dividing by the 
uncorrected G star spectrum, and multiplying the result 
with the esdM spectrum. These corrections were gener- 
ally adequate, but noticeable residuals are seen around 
the O2 A-band in the spectra of SSSPM 0500-5406 and 
APMPM 0559-2903, for which G star telluric cafibra- 
tors were observed at a larger differential airmass. These 
residuals principally affect the blue wing of the K I dou- 
blet at 7665/7699 A. 

Reduced spectra for the four esdMs observed are shown 
in Figure|21 Each spectrum shows characteristic spectral 
traits of late- type esdMs, including weak TiO absorp- 
tion bands at 6400, 6700 and 7150 A (and possibly weak 
TiO absorption at 7800 A); strong CaH bands at 6400 
and 6900 A; and prominent line absorption from K I 
(7665/7699 A doublet) and Na I (8183/8195 A doublet). 
The K I doublet lines show substantial broadening in the 

* IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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spectrum of LEHPM 2-59, reminiscent of the K I broad- 
ening observed in solar metallicity L dwarfs. The 6708 A 
Li I hne is not detected in any of the spectra; this species 
has almost certainly been depleted by fusion reactions in 
the cores of all of these objects. There is a dip in the 
spectra of SSSPM 0500-5406, APMPM 0559-2903 and 
LEHPM 2-59 at 7400 A which we attribute to the 1-0 
B^'n-X'^S" VO band, a feature commonly observed in 
solar-metallicity M dwarf and early L dwarf spectra but 
identified here for the first time in esdM spectra. This de- 
tection is robust despite poor telluric O2 correction in the 
spectra of SSSPM 0500-5406 and APMPM 0559-2903, as 
this calibration error affects only the 7600-7700 A win- 
dow, while a dip in the spectrum is seen clearly at 7400- 
7500 A (se e also 5). Using the Kurucz Atomic Line 
Database^ ijKurucz fc Be ll 1995), we have also identified 
a forest of Ca I lines between 6100-6170 and 6440-6500 
A, as well as more prominent lines at 6573 and 7326 A. 
All of the atomic features identified in these spectra are 
listed in Table El molecular line identifications can be 
found in Kirkoatrick. Hcnrv, & McCarthv ( 1 991') . 

Rb I lines at 7800 and 7948 A are also present in the 
spectra of APMPM 0559-2903 and LEHPM 2-59, and the 
7948 A line is weakly present in the spectra of LP 589-7 
and SSSPM 0500-5406. These lines are particularly in- 
teresting, as they have been previously observed o nly in 
the lo w resolution spectra of L dwarfs l|Ki rkDatric k et al.l 
119991) a nd, more recently, of the esdM6.5 LSR 0822-1700 
(jLepinc . Shara, & Rich 20 0^ and the L subdwarf LSR 
1610-0040 ()Lepine. Rich, fc Shar a 2003b, however, see 
Gushing fc Vacca 2006; Reiners fc Basri 2006). On the 
other hand, Rb I is not detected in low resolution spectra 
of the latest-type M dwarfs or sdMs l'Kirk£gt^ck_gj^n 
[1299; Lepine. Shara. & Rich 2003; Schol z et alJl2004aD . 
This atomic species is largely ionized in the photospheres 
of M dwarfs (T„„ ~ 2350 K; Lodders 1999, 2002), al- 
though a significant fraction of Rb I gas is probably neu- 
tral at higher temperatures in the higher pressure photo- 
spheres of subdwarfs. The appearance of both the 7400 
A VO band and the 7800/7948 A Rb I doublet lines could 
be considered a defining trait of ultracool esdMs. 

Finally, we note that Ha emission, a feature arising 
from magnetic activity that is nearl y alway s present 
in the spectra of late-tvpe M dwarf s l|Gizis et al. 2000; 
iHawlev et al.ll2002t lv\^st et al.l l2004l'). is not detected in 
any of these spectra (including LP 589-7; Gizis & Reid 
1999), consistent with weak or absent chromospheres. 

4. ANALYSIS 

4.1. Atomic Line Strengths and Radial Velocities 

Equivalent widths (EWs) for the alkali lines and the 
two most prominent Ga I lines were measured using the 
IRAF SPLOT routine. Uncertainties were determined as 
the scatter of values from multiple measurements, while 
upper limits for non-detected lines (Rb I in LP 589-7 and 
SSSPM 0500-5406) were estimated from the mean EWs 
of local noise features. These values are given in Tabled 
Note that EWs given for the K I and Na I doublets are 
the combined absorption from both lines in each feature. 
Nearly all of these lines are seen to strengthen from LP 

^ Obtained through the online database search form 
created by C. Heise and maintained by P. Smith; s ee 
http : //cf a-ww¥ . harvard . edu/amdata/ ampdata/kurucz23/ sekur . html 



589-7 to LEHPM 2-59, with the possible exception of the 
7326 A Ca I line, a transition that has the highest lower 
energy state in Tabled (4.6 eV; Kurucz & Peytremann 
1988). The strengthening of atomic lines is observed in 
the spectra of solar-metallicity M and L dwarfs as one 
progresses to later spectral types, and is attributed to 
lower photospheric temperatures and the corresponding 
increase in column abundances of these neutral species. 
The evolution of the line strengths in the esdMs observed 
here also suggests a trend of cooler T^f / with later esdM 
type. 

The six alkali lines and two strongest Ga I lines were 
used to measure radial velocities (Vrad) for the observed 
sources. Line centers were determined by gaussian fits 
to the line cores, and velocity shifts were determined rel- 
ative to the vacuum wavelengths listed in the Kurucz 
database (Table 131 . The mean and standard deviations 
of these shifts are listed in Table 01 and include a sys- 
tematic uncertainty of 3 km s~^ based on the mean accu- 
racy of the dispersion solutions. Our values for LP 589-7, 
SSSPM 0500-5406 a nd APMPM 5 59-29 03 are in agree- 
ment with those of iGizis fc ReidI l)1999j) : iLodieu et alJ 
(|2005f) : and lSchweitzer et alJ l|199l . respectively. With 
Vrad = 79±9 km s~\ LEHPM 2-59 also appears to have 
kinematics consistent with a halo subdwarf. 

4.2. Spectral Classification 

The optical spectra were classified using the scheme 
of G97 as updated by LRS03. This scheme is based 
on the relative strengths of the 7100 A TiO band 
and the 6400/6900 A GaH bands, as sampled by the 
four indices GaHl, G a H2, G aH3 and Ti05 defined in 
iReid. Hawlev. fc GizisI l) 19951) : as well as the redness 
of the pseudo-continuum in the 6500-8200 A spectral 
region as sampled by the Golor-M index of LRS03. 
These indices were measured for each of the spectra af- 
ter shifting them to their rest frame velocities. Fig- 
ure 0] compares the combined GaH2-|-GaH3 indices to 
Ti05 for these sources, in addition to measurement s 
for late-type dwarfs from Hawlev. Gizis. fc Reid (1996*); 
G97; Gizis fc R eid (1997); Reid & Gizis (2005); LRS03; 
iLepine. Shara. fc Rich (2004); and Scholz et al. (2004a). 
This index-index diagram, based on work by G97, has 
been used by several groups to segregate M dwarf spectra 
by metallicity class. We update the divisions originally 
set forth by G97 for the combined GaH2-|-GaH3 indices: 
sdM : (GaH2 + GaH3) < 1.31(Ti05)3 - 

2.37(Ti05)2 + 2.66(Ti05) -0.20 (1) 

esdM : (GaH2 + GaH3) < 3.54(Ti05)3 - 

5.94(Ti05)2 + 5.18(Ti05) - 1.03. (2) 

We stress that, as in G97, these divisions are somewhat 
arbitrary as sources spanning Figure ^represent a broad 
continuum of metallicities. Figure 0] nevertheless demon- 
strates that SSSPM 0500-5409, APMPM 0559-2903 and 
LEHPM 2-59 are aU bona-fide esdMs lying at the tail 
end of the esdM distribution. LP 589-7 appears to be a 
borderline sd/esd object; however, for the remainder of 
this paper we will treat it as an esdM. 

Numerical subtypes for the four observed sources were 
determined using the relations 

SpTesdM = 7.91(GaH2)2 - 20.63(GaH2) + 10.71 (3) 
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SpTg^dAf = -13.47(CaH3) + 11.50 (4) 

(i.e., Eqns. 2 and 8 from G97), and 

SpTesdM = 22.0 ln(Color - M) - 4.1 (5) 

(i.e., Eqn. 146 in LRS03), where SpTe^dM(esdMO) = 0, 
SpTesdM(esdM5) = 5, etc. Technically, these relations 
are defined only up to subtype esdM6 (G97) , but we fol- 
low current practic e (Lcoinc. Shara. & Rich 200?, '2004'; 

IScholz et alJl2004at .Lodieu et al 200 51 in extrapolating 

these relations to later numerical types (however, see 
§6.2). Index values, index subtypes and averages of the 
numerical subtypes (rounded off to the nearest 0.5 sub- 
type) are listed in Table |S1 

The derived subtypes for LP 589-7 and APMPM 
0559-2903 are consistent with previous determinations. 
We derive an esdM6.5 classification for SSSPM 0500- 
5406, 0.5 sub t ypes later than the classification given by 
iLodieu et al.l (j2005) who measured CaH2 and CaH3 in- 
dices 0.05 larger than our values. However, this is gen- 
erally consistent with the typical uncertainties in index 
measurements and the corresponding uncertainties in de- 
rived spectral types (±0.5 subtype). For LEHPM 2-59 we 
derive a classification of esdM8, a full subtype later than 
APMPM 0559-2903, until now the latest esdM known. 
This result is consistent with the observed spectral trends 
noted above, particularly the deeper near infrared H2O 
bands and stronger optical atomic lines in the spectrum 
of LEHPM 2-59. We conclude that LEHPM 2-59 is the 
latest-type esdM found to date. 

4.3. Distance and Space Velocity Estimates 

There are currently few late- type esdMs with measured 
parallaxes; only two sources esdM5 and later (the esdM5 
LHS 3061 and the esdM6 LHS 1742a; Monet et al. 1992) 
have known distances. Hence, only a rough distance es- 
timate {dest) for LEHPM 2-59 is possible. For this, we 
employed the spectral type/absolute magnitude relations 
of LRS03 (Eqns. 24 and 25), which can be rewritten as 

Mfl = 10.44 -h0.49(SpT,,,,,) (6) 

= 8.06 + 0.37(SpT,,rfM), (7) 

where R is from USNO-Bl.O iMonet et al.ll2003(l and Ks 
is from 2MASS. iReid fc GizisI l)2005|) derived a similar 
relation for Mk^ versus spectral type for esdMs, 

Mk. =8.73 + 0.31 (SpT,,,m)- (8) 

Assuming SpTe^dM = 8, i? = 18.82 and = 14.76 (Ta- 
ble [TJ, these three re lations yield d i stance s of 79, 55 and 
51 pc, respectively. iGizis fc ReidI ((X999') have also de- 
rived absolute magnitude/color relations for esdMs; in 
particular, 

M/ = 5.96 + 4.29(i?-/), (9) 

where R and / are photographic magnitudes. Again, 
adopting photometry from USNO-Bl.O [R = 19.92, / 
= 17.24) yields dest = 79 pc, in agreement with the 
Mij/spectral type relation in Eqn. 6. We therefore adopt 
an average dest — 66±15 pc for this source, but stress 
that such estimates are highly uncertain given the cur- 
rent paucity of parallaxes for late-type subdwarfs. 

" Note that the LRS03 Color-M/spectral type relation uses the 
natural logarithm of the Color-M index, and not the base-10 loga- 
rithm as suggested by their Eqn. 14. 



Using this estimated distance, and the measured 
proper motion and radial velocity of LEHPM 2-59, we 
computed UVW space velocity components relative to 
the Local Standard of Rest (LSR). We derive [U,V,W\ 
« [135,-180,-80] km s~^, assuming an LSR sol ar motion 
oUUq.Vq.Wq] = [10,5,7] km s'^ (Do hTicn Binnevl 
11998ft . These velocities lie well outside of the velocity 
dispersion sphere of local disk M dwarfs ([tT[/,tTy ,(Tvf] ~ 
[40,28,19] km s^^ centered at [-13,-23,-7] km s^^ Haw- 
ley, Gizis & Reid 1996); and the large negative V veloc- 
ity, ranging over -210 to -145 km s^^ for our distance 
estimates, is consistent with motion independent of the 
Galactic disk . Assu ming an LSR rest frame velocity of 
-220 km s"^ l|Kerr fc Lvndcn-BcU 198^), the total space 
velocity of LEHPM 2-59 in the Galactic potential is 135- 
190 km s"^. This is well below the Galactic escape veloc- 
ity in the vicinitv of the Sun ijCarnev. Latham. &: LairdI 
ll988HLeonard fc Tremaindll990|) . Interestinglv. the or- 
bit of this object is highly elliptical, and passes well inside 
the Galactic bulge on its closest approach to the Galactic 
center. This is not atypical for halo stars currently in the 
vicinity of the Sun (S. Lepine, 2006, private communica- 
tion). 

5. SPECTRAL MODEL FITS 

To further gauge the physical properties of these late- 
type esdMs, we compared our optical and near infrared 
spectra t o subsolar metallicity theoret i cal sp ectral mod- 
els from iHauschildt. Allard fc BaronI l)1999l NextGen) 
and lAllard et al.ll|200lL AMES Gondl. These models are 
based on the Phoenix code fHauschildt. Baron & AUar^ 
119971 and references therein) , employ self-consistent tem- 
perature/pressure profiles, and assume local thermody- 
namic equilibrium. Both sets of models use the opac- 
ity sampling method; a full account of the chemical 
sp ecies and opacities used by these models is provided 
in lAUard et alJ l)200lD and references therein. Two ma- 
jor changes in the Cond models are the removal of con- 
densed species from the upper atmosphere and the use of 
updated TiO and H2O opacities. The NextGen models 
and their antecede nts (Allard 1990; Allard & HauschildJ 
1995; Allard et al. 1997) have been used extensively for 
fittin g M subdwarf spectra [Gizis. 1997; .Schweitzer ct all 
799^ jPawson De RobertiFi 120001 iLeggett et a,1.ll200nt 
Lepine. Shara. & Richll200l . 

We sampled grids of NextGen and Cond models span- 
ning temperatures of 2600 < Tey/< 3600 K in steps of 
100 K, metallicities of -3.0 < [M/H] < 0.0 dex in steps 
of 0.5 dex, and surface gravities of log 5 = 5.0 and 5.5 
cm s~^. Comparisons were made separately to the op- 
tical and near infrared data. For the optical fits, both 
empirical and model spectra were normalized at 8100 A 
as in Figure 13 and the observed data were shifted to 
their rest frame velocities. For the near infrared fits, 
data and models were scaled to 1.2 /im. Model spectra 
were also deconvolved to the resolution of the observed 
data (A/AA = 1800 and 150 for the optical and near in- 
frared, respectively) using a Gaussian kernel. For each 
spectrum/model pairing, the root mean square (RMS) 
deviation was computed over the spectral ranges 6050- 
7500 and 7650-8250 A in the optical (to avoid regions of 
poor telluric O2 correction and the 7400 A VO band; see 
below); and 0.75-1.31, 1.45-1.75 and 2.0-2.4 ^m in the 
near infrared (to avoid telluric H2O regions). The nor- 
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malization of the model spectra in each fit was allowed 
to vary by a factor ranging from 0.75 to 1.25 (in steps 
of 0.05) to allow for continuum offsets, and the normal- 
ization with the minimum RMS was retained. For each 
model set and surface gravity, the Teff and [M/H] com- 
bination with the overall minimum RMS was deemed to 
be the best fit for a particular spectrum. These best fit 
values are listed in Table |S1 

Figure El and [HI compares the best fit NextGen and 
Cond, log 5 = 5.0 models for the spectrum of LEHPM 
2-59 in the optical and near infrared, respectively. Over- 
all, these fits are reasonably good, with the broad spec- 
tral shape and several individual features matching quite 
well. For the optical data, the log 17 = 5.0 NextGen mod- 
els consistently provided the best fits to the data. This 
appears to be largely due to better fits to the 6700 A CaH 
band which is only slightly overestimated in the NextGen 
models but greatly underestimated in the Cond mod- 
els. The 7150 A TiO band is too strong in both model 
sets; higher Teff and/or lower [M/H] models weaken this 
band but show very poor agreement with the rest of the 
optical spectrum. TiO features have long been identi- 
fied as a p roblem, particularlv in the NextGen models 
ijAUard. H auschildt fc Schwenkell2000() . Atomic lines are 
also generally too strong in the best fit models. Indeed, 
there are many metal lines apparent in the Cond mod- 
els for which no analogs are seen in the empirical data. 
The Na I, Rb I and prominent Ca I lines are also too 
deep, although the K I doublet is reproduced quite well, 
both in depth and breadth, for both sets of models. The 
apparent mismatch over 7300-7500 A is due to the 7400 
A VO band, which is not included i n the opacity set o f 
either the NextGen or Cond models l|Allard et al.ll20niD . 
We found slight systematic effects in the derived param- 
eters between different model sets and surface gravities. 
The best fit Cond models yielded slightly lower Te//s 
(by < 100 K) and higher [M/H] values (by < 0.5 dex) as 
compared to the best fit NextGen models; higher gravity 
model fits typically gave higher Te//s (by < 100 K) and 
[M/H] values (by < 0.5 dex). 

In the near infrared, spectral model fits again provided 
reasonably good matches to the broad spectral energy 
distribution, with slight discrepancies in the 0.99 /im FeH 
band, the 1.15-1.30 /im continuum (also dominated by 
FeH absorption in late-type subdwarfs; Burgasser et al. 
2004, Gushing et al. 2006; Burgasser et al. in prep.) and 
the 1.4 /im H2O band. In this wavelength regime, Cond 
models typically provided the best fits, largely due to 
better agreement in the 0.8-1.1 fim region. Systematic 
shifts in the derived parameters between the NextGen 
and Cond model sets were less pronounced at these wave- 
lengths, but higher surface gravity models consistently 
gave larger Te//s (by 100-200 K) and shghtly lower 
[M/H] values (by < 0.5 dex). There was also some degen- 
eracy in the best fit values, with hotter, solar metallic- 
ity models providing fairly good matches to the data for 
A > 1.1/im. This is not surprising, given that the blue, 
relatively featureless near infrared spectra of esdMs are 
not unlike those of hotter M and K dwarf stars. 

Derived parameters between the optical and near in- 
frared fits for a given object exhibit clear systematic dif- 
ferences. TeffS derived from the optical data are gener- 
ally 100-200 K higher than those derived from the near 
infrared data; metallicities are consistently 0.5-1.0 dex 



higher. Examination of the models indicates that lower 
metallicities are required to match the H- and K-hand 
suppression in the near infrared data, but can result in 
optical CaH bands that are too deep. Similar effects 
are observed with lower T^ffS. Since we have no way 
of independently validating fits in either spectral region, 
we must treat these systematic differences as a source of 
uncertainty in the model fits, emphasizing that absolute 
Te/ys and [M/H] values derived from spectral model fits 
should in general be treated with some caution. 

Turning to the derived parameters of the objects in 
our sample, we find that, overall, later esdM types cor- 
respond to cooler Te//s, consistent with expectations. 
The one notable exception is the esdM6.5 LSR 0822- 
1700, for which we derive a temperature 100 K cooler 
than the esdM7 APMPM 0559-2903. This shght dif- 
ference may not be significant; indeed, our Teff deter- 
mination based on near infrared data is 200 K cooler 
than that derived bv lLepine. Shara. fc RichI l)2004D from 
optical data. However, metallicity effects may be at 
play. LSR 0822-1700 appears to have a consistently lower 
[M/H] than APMPM 0559-2 903 based on both optical 
ijLepine. Shara. fc RichI 120041) and near infrared analy- 
sis. In any case, LEHPM 2-59 is the lowest temperature 
object in the group, with Teff ^ 2800-3000 K. Evolu- 
tionary models by Baraffc et al. (1997) predict a mass 
of ^0.09 M0 for these parameters assuming an age of 
10 Gyr, only ~0.01 M© above the hydrogen burning min- 
imum mass for these metallicities. 

6. DISCUSSION 

6.1. The Temperature Scales of M Dwarfs and Extreme 

Subdwarfs 

The temperatures^ d erived for late-type esdMs in this 
and other studies (j Schwcitz er et aLlll999l: iLeeeett et alJ 
120001 iLepine. Shara. fc Rich<]2004D are relatively warm 
compared to solar-metallicity dwarfs with comparable 
subtypes, which typically have T^jf ^ 2200-2800 K 
(KirkDatrick et al. 1993; Lcgge ^^AU^^^ fc HauschildJ 
lT998t il'l^^t era].lf200(Tt iGolimowski e Tal .1 1200411 . Fig- 
ure[7|compares Te// determinations for M dwarfs and es- 
dMs based on spectral model fits to NextGen models and 
their antecedents. These data are segregated by optical^ 
^Cirkpatrick et al. 1993; Gizis 1997; Schweitzer et a^ 
'1999; Lcpine, Shara, fc Rich 2004) and near infrared 

( Kirk patrick et aL.1993. ..Legeett et al 1996. ,2000- ,20011 

Daws on fc De R,obertisl 120001) analvses. There is a sig- 
nificant downturn in the dwarf Tef / scale around type 
M5-M6, dropping from -3000 K to ~2200 K by type 
M8. The esdMs, on the other hand, show a more grad- 
ual decrease in Teff with spectral type over the same 
range. Linear fits to optical and near infrared Teff de- 
terminations over the range esdMl-esdM8 yield parallel 
relations: 

Teff = 3750 - 93xSpT^,rf^,f K (optical) (10) 

Teff = 3620 - 96xSpT^,rfA,j K (near infrared) (11) 

with dispersions of 30 and 90 K, respectively. The tem- 
peratures of the two latest esdMs, APMPM 0559-2903 

^ We did not include T^ff determinations for M dwarfs later 
than M6 in the |KirkDatrick ct al. (1993) and G97 studies due to 
the poor quality fits to ultracool dwarf spectra existing models 
provided at that time. 
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and LEHPM 2-59, are 600-800 K hotter than those of 
equivalently typed solar- metalHcity dwarfs. 

Why would the effective temperatures of ultracool es- 
dMs differ so greatly from those of ultracool solar metal- 
licity dwarfs? It is well known that lower metallicity 
implies lower overall atmospheric opacity, a consequence 
of reduced molecular and H~ abundances. Hence, for a 
given mass, the observed photosphere lies at deeper, hot- 
ter and higher pressure layers in more metal-poor stars. 
However, this argument serves only to explain (possi- 
ble) differences in Te///mass relations between metal- 
hcity classes. Te// /spectral type relations also hinge 
on how the spectral types themselves are derived. In 
fact, the classification of M dwarfs and esdMs are de- 
termined by separate relations in the G97 and LRS03 
schemes, as necessitated by the divergence of spectral 
properties between these metallicity classes. Hence, dif- 
ferences in the Tg// /spectral type relations between M 
dwarfs and esdMs may simply be an artifact of the un- 
derlying classification scheme. At first glance, this may 
seem unimportant, for as long as the correct scale is used 
for a given metallicity class, one should derive the cor- 
rect Teff- However, metallicity is not a discrete quantity 
like metallicity class, and systematic deviations can occur 
for individual objects that are more or less metal-poor 
than the mean of their assigned class. A more global 
Te// /metallicity /spectral type relation is required. 

6.2. Future Prospects for Ultracool Subdwarf 
Classification 

With the recent compilation of several red optical 
proper mo tion surve ys and the initiation of near infrared 
programs l|Deacon. Hamblv & Cooke 2005, Looper et al., 
in prep.), it is quite likely that LEHPM 2-59's status as 
the coolest esdM known will not be long-lived. Indeed, 
the identification of a few L subdwarfs (sdL) to date sug- 
gests that esdL discoveries may not not f ar off (however , 
see below). It has been show n by Burgass er et alJ l)2003af) 
and lLepine. Rich, fc Sharal (|2003b.') that these ultracool 
metal-poor stars and brown dwarfs cannot be adequately 
classified by existing schemes largely due to the disap- 
pearance of the 7150 A TiO band. The lack of flux from 
these cool objects also argues against classification in the 
6300-7200 A region. 

So where do we proceed? An obvious option is to de- 
fine new subdwarf classification schemes at longer wave- 
lengths, including the red optical (cf. late-M and L 
dwarfs) and near infrared (cf. T dwarfs) spectral regions. 
Absorption bands from the metal hydrides FeH and CrH 
(and perhaps TiH; Burrows et al. 2005) are prominent in 
ultracool dwarf and subdwarf spectra for A > 8000 A, and 
are likely to be present in even lower metallicity analogs. 
Contrasting these features with longer wavelength TiO 
(8200 and 8400 A) and VO bands (7900, 9500 and 10500 
A), alkali lines (Na I, K I, Kb I and Cs I) or pseudo- 
continuum slope should prove to be effective diagnos- 
tics for segregating metallicity and temperature classes. 
At longer wavelengths, subdwarf classification becomes 
more difficult due to the weakness of absorption features 
and the strong suppression of H- and K-hand flux by H2 
absorption (which also wipes out the 2.3 fim CO band). 
Nevertheless, ratios comparing the 1.4 /im H2O band and 
the near infrared spectral slope might provide some dis- 



crimination (Burgasser et al. in prep.), while metal line 
features (including Al I; see Cushing & Vacca 2006) and 
FeH bands in the I-I.35 /im region could be used for 
higher resolution studies. 

As later extreme subdwarfs are identified, what will 
define the termination of the esdM class and the begin- 
ning of the esdLs? Solar-metallicity L dwarfs are dis- 
tinguished by waning metal oxide bands, strengthening 
metal hydride and alkali lines, steep red optical slopes 
and red near infrared colors. Yet esdMs already exhibit 
weak metal oxides and strong metal hydrides, and sev- 
eral neutral alkali species (including Rb I and Cs I) ; while 
near infrared color s will never become r ed due to strong 
absorption by H2 ijSaumon et alJlT99^ . Spectral mod- 
els cannot provide clear guidance on this transition, as 
it remains unclear as to whether condensate formation, 
a key aspect of the M/L dwarf transition, plays a sig- 
nificant ro le in Iqw-temp erature metal-poor atmospheres 
(|Burgasscr et al.ll2003af) . Clearly, cooler extreme subd- 
warf discoveries must guide the definition and our under- 
standing of this transition. 

One further complication in this issue is whether halo 
esdLs even exist in our Galaxy today. As metal-poor 
halo objects about the substellar mass limit evolve over 
'-~^I0 Gyr, wide gaps in the luminosities and effective 
temperatures of this population develop. Low mass 
stars attain a steady-state T^ff > 2500-3000 K, while 
most brown dwarfs c ool to Te// < 1000 K in this time 
(jBurrows et al.l 120011 cf. Figure 3 in Burgasser et al. 
2003a). Only a narrow range of low mass star/brown 
dwarf transition objects are expected to encompass in- 
termediate L-type Te//s between these limits, and could 
therefore be quite rare. Future surveys for low luminosity 
metal-poor objects will hopefully probe this L dwarf gap, 
providing, if little data for esdLs, a unique constraint for 
brown dwarf evolutionary theories. 

Finally, one area in which subdwarf classification can 
be immediately improved is the identification of spe- 
cific standard stars for temperature and metallicity sub- 
types. A framework of standard stars is a fundamental 
tenet of the MK Process (Morgan. Kecnan & Kellman 
'1943' 'Morgan & Keenan 1973; Keenan fc McNeill Il97a 
iCorballv. Grav fc Garrison .199^ the most widely 
adopted method of stellar classification. Standard 
stars provide the basis of current classification schemes 
for so lar met allicity M ( Kirkoatrick. Hcnrv. fc McCarthy 
[l99l|l). L (Kirkpatric k et alJ a999) and T dwarfs 
(jBurgasser et al. .2006) . Allowing specific stars to de- 
fine a classification scheme provides a level of consistency 
that is not generally present in pure index schemes, while 
retaining independence from constantly evolving theo- 
retical interpretations. There are now several dozen M 
subdwarfs of types sdM5/esdM5 and earlier from which 
appropriate standards can be chosen, and we anticipate 
that current red optical and near infrared proper motion 
surveys will soon fill in the remainder of the M subdwarf 
and extreme subdwarf sequences. Now is an opportune 
time to consider the construction of a more robust clas- 
sification scheme for late-type subdwarfs. 

7. SUMMARY 

We have identified LEHPM 2-59 as the coolest esdM 
identified to date. Near infrared and optical spectroscopy 
show features indicative of a late-type esdM, including 
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strong alkali lines, VO and H2O absorption and a blue 
near infrared spectral slope. We derive an optical spec- 
tral type of esdMB on the G97 and LRS03 schemes, and 
spectral model fits to both optical and near infrared data 
indicate Tg// = 2800-3000 K and -1.5 < [M/H] < -2.0 for 
this source. Its kinematics confirm it as a halo star. We 
show that the temperatures of this and other late-type 
esdMs are significantly higher than equivalently classi- 
fied solar-metallicity M dwarfs, a divergence that could 
lead to systematic errors in derived parameters for ob- 
jects with intermediate metallicities. Finally, we have 
touched on the future prospects of M subdwarf classifi- 
cation, and methods by which it may bo improved and 
extended to later subtypes. With the current high dis- 
covery rate of low-temperature, ultracool halo subdwarfs, 
it is likely that these issues will require further scrutiny 
in the near term. 
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Fig. 1. — Hj versus (Reso ~ J) diagram for the LEHPM catalog IPokornv et all2004l) . Note the three clusters of sources corresponding 
to (from center to lower left) main sequence dwarfs, halo subdwarfs and white dwarfs. Our selection criteria for ultracool dwarf and 
subdwarf candidates are indicated by dashed lines. Previously classified sources within this region are noted by open circles; LEHPM 2-59 
is indicated by an open square. [THIS FIGURE IS INCLUDED IN JPEG FORMAT] 




1.0 1.5 2.0 

Wavelength (/i.m) 



Fig. 2.— SpeX prism spectra of the esdMs LP 589-7 (esdMS), LSR 0822+1700 {esdM6.5), APMPM 0559-2903 (esdM7) and LEHPM 
2-59 (esdM8). Spectra are normalized at their flux peaks and offset by constants (dotted lines). Prominent features in the NIR spectra of 
cool esdMs discernible at the resolution of these prism data (A/ A A 150) are labeled. Regions of high telluric opacity are indicated by ffi 
symbols. 
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Fig. 3.— Optical spectra of the late esdMs LP 589-7 (esdM5), SSSPM 0500-5406 (esdM6.5), APMPM 0559-2903 (esdM7) and LEHPM 
2-59 (esdM8) from top to bottom. All data a re shifted to their rest frame velocities, n ormalized at 8100 A and of fset by a constant (do tted 
lines). Identified molecular features are from lKirkpatrick. Henry. &: McCarthy! 11991i) : atomic features are from lKurucz fc Bell tT99l see 
Table 3). Residual noise from the telluric O2 A-band is indicated by © symbols. 
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Fig. 4. — Spectral indices CaH2+CaH3 versus Ti05 for dwarfs (points), subdwarfs (open squares) and extreme subdwarfs (filled triangles) 
fromlHawle v. Gizis. & Reid (1993); G97; pizis & Reid (1997); Reid et al. (2002); LRS03; Lcpine. Shara. & Rich (2004); and Scholz et a]] 
J200la"). Data from this paper are encircled and labeled, while values for the esdM6.5 LSR 0822-1-1700 and the sdM9.5 SSSPM 1013-1356 
from Lcpine. Shara. & Rich (2004) and Scholz et al. (,2004al) . respectively, are also labeled. Dashed lines delineate boundaries between 
dwarfs, sdMs and esdMs as defined by Egns. H^andl^ 
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Fig. 5. — Comparison of best fit NextGen (top) and AMES Cond (bottom) \ogg = 5.0 cm s ^ models to the observed red optical 
spectrum of LEHPM 2-59 (black). LEHPM 2-59 spectra are normalized at 8f00 A; model spectra are scaled to their best fit normalization. 
The wavelength ranges for which spectral data and models were compared are indicated by the hatched areas. 
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Fig. 7. — Tef fS for field M dwarfs (circles) and esdMs (squares) based on spectral model fits. Data are from'Kirkoatrick ct al.' ^1991 
SpT < M6l: ILeggett et alJ tl996. 2000, 200ll): IGizT3 fl997,. SpT < M6l : ISchweitzer et al.l tT99l ): iLepine. Shar a. & Rich (2004); and this 
paper (oversize symbols). Tf,//s derived from fits to optical data are indicated by open symbols, those from near infrared data by solid 
symbols. The solid line delineates the M dwarf Teff scale from lR.eid fc Hawlev 1.2000.') i the dashed lines delineate linear fits for esdM T^ffS 
based on optical (top) and near infrared (bottom) data (Egns. [THl and lTTl . 
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TABLE 1 
SpeX Observing Log 



Source 



UT Date tint (s) Airmass Flux Cal. SpT Ref/ 



LP 589-7 
LEHPM 2-59 
APMPM 0559-2903 
LSR 0822-1700 



01'' 57*" 27^92 -1-01° 16'" 43'.'3 14.50±0.03 2004 Sep 5 1080 1.13 HD 13936 AO V 1 

04h 09"94 -22° 45'" 08'.'4 15.52±0.05 2004 Sep 9 720 1.43 HD 32855 AO V 2 

05'' 58" 58"91 -29° 03"^ 26'.'7 14.89±0.04 2005 Dec 31 1440 1.53 HD 41473 AO V 3 

08'' 22™ 33^69 -17° 00"' 19'.'9 15.87±0.08 2004 Mar 10 1080 1.01 HD 58383 AO V 4 



References. — (ll KJizia heij ITSSai ; f2> IPokornv et aUHTOtl ^ f 31 mdiweitzer et alJimiiill ^ f4l lLepine. ijkara. fc U-lcU IW!^ 

^,12000 Coordinates from 2MASS. 
^ J magnitudes from 2MASS. 
^Discovery reference for esdM souree. 



TABLE 2 
LDSS-3 Observing Log 



Source 










5^ 






UT Date 


tint (s) 


Airmass 


Tell. Cal. 


SpT 


Ref.<: 


LP 589-7 


01'' 


57" 


27'.'92 


+01° 


16™ 


43"3 


17.39 


2005 Dec 4 


2100 


1.16 


HD 603 


G2 V 


1 


LEHPM 2-59 


04'' 


52™ 


09'.' 94 


-22° 


45™ 


08'.' 4 


18.72 


2005 Dec 4 


1800 


1.04 


HD 31527 


G2 V 


2 


SSSPM 0500-5406 


05'' 


00™ 


15'.' 77 


-54° 


06™ 


27'.' 3 


17.42 


2005 Dec 4 


600 


1.50 


HD 33967 


G2/3 V 


3 


APMPM 0559-2903 


05'' 


58™ 


58^91 


-29° 


03™ 


26'.' 7 


18.08 


2005 Dec 4 


1200 


1.06 


HD 33967 


G2/3 V 


4 



REFERENCES. — fl^ lUizis j^. ^^7^111??^ : f2) IPokornv eT^ fS^ ILoJieu e^"^ nHH^ ; f 4) IblcWeUzer aUIT?^ . 



^,12000 Coordinates from 2MASS. 

^Photographic R magnitudes from the SupcrCosmos Sky Survey. 
'^Discovery reference for esdM source. 



TABLE 3 

Atomic Lines Identified in esdM Optical Spectra 
OVER AA 6100-8300 A 



Wavelength Element 


^lower 


E upper Ptcf. 


(A) 


(eV) 


(eV) 



6102.723 


Cal 


1, 


,879467 


3.910663 


1 


6122.217 


Cal 


1, 


.885935 


3.910663 


1 


6162.173 


Cal 


1, 


.899063 


3.910663 


1 


6166.439'' 


Cal 


2, 


,521433 


4.531641 


1 


6169.042" 


Cal 


2, 


,523157 


4.532517 


1 


6169.563" 


Cal 


2, 


,525852 


4.535043 


1 


6439.075 


Cal 


2, 


,525852 


4.450947 


1 


6449.808 


Cal 


2, 


,521433 


4.443325 


1 


6462.567 


Cal 


2, 


,523157 


4.441254 


1 


6471.662 


Cal 


2, 


,525852 


4.441254 


1 


6493.781" 


Cal 


2, 


,521433 


4.430310 


1 


6499.650" 


Cal 


2, 


,523157 


4.430310 


1 


6572.779 


Cal 


0, 


,000000 


1.885935 


1 


6717.681 


Cal 


2, 


,709192 


4.554447 


1 


7148.150 


Cal 


2, 


,709192 


4.443325 


2 


7202.200 


Cal 


2, 


,709192 


4.430310 


2 


7326.145 


Cal 


2, 


,932710 


4.624710 


2 


7664.911 


KI 


0, 


,000000 


1.617220 


1 


7698.974 


KI 


0, 


,000000 


1.610064 


1 


7800.259 


Rb I 


0, 


,000000 


1.589158 


3 


7947.597 


Rb I 


0, 


,000000 


1.559697 


3 


8183.255 


Nal 


2, 


,102439 


3.617221 


4 


8194.790" 


Nal 


2, 


,104571 


3.617221 


4 


8194.824" 


Nal 


2, 


,104571 


3.617215 


4 



References. — (1) 'Wies e^ S mith, fc GlennorJ 

IKurucz fc Pevtremanni (il975il . 
^Blend. 
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TABLE 4 

Alkali Atomic Line Equivalent Widths (A). 



Source 


Cal 


Cal 


KI 


Rbl 


Rbl 


Nal 


Vrad 




(6573 A) 


(7326 A) 


(7665/7699 A) 


(7800 A) 


(7948 A) 


(8183/8195 A) 


(km s-1) 


LP 589-7 


1.17±0.06 


0.76±0.08 


9±2 


<0.10 


0.28±0.03 


4.54±0.14 


-82±9 


SSSPM 0500-5406 


2.03±0.15 


0.99±0.16 


26±2 


<0.15 


0.5±0.2 


6.5±0.4 


216±9 


APMPM 0559-2903 


2.0±0.2 


1.1±0.2 


27.7±1.9 


0.34±0.10 


0.62±0.11 


6.74±0.14 


181±5 


LEIIPM 2-59 


2,]7±n.l3 


0.93±n.l2 


39.7±1.2 


o.60±n.n7 


0.68±n.05 


7.2±0.3 


79±9 



TABLE 5 
Spectral Ratios and Classification. 



Source 



CaHl 



CaH2^ 



CaH3'^ 



Ti05 Color-M'» 



SpT 



LP 589-7 0.472 0.325 (4.8) 0.493 (4.9) 

SSSPM 0500-5406 0.310 0.220 (6.5) 0.331 (7.0) 

APMPM 0559-2903 0.342 0.217 (6.6) 0.331 (7.0) 

LEHPM 2-59 0.267 0.175 (7.3) 0.265 (7.9) 



0.649 1.49 (4.7) sd/esdM5 

0.755 1.64 (6.8) esdM6.5 

0.604 1.67 (7.1) esdM7 

0.656 1.79 (8.7) esdMS 



''Numerical esdM subtype in parentheses based on Eqns. 1-3. 



TABLE 6 
Spectral Model Fits.* 



Source 



SpT 



(K) 



NextGen 
5.0*' 5.5 
P7HJ 



(dex) 



(K) 



W7W 

(dex) 



(K) 



AMES Cond 
5.0 5.5 



(dex) 



(K) 



W7W 

(dex) 



Optical 



LP 589-7 


sd/csdM5 


3300 


-1.0 


3300 


-1.0 


»200 


-1.0 


3300 


-1.0 


SSSPM 0500-5406 


csdM6.5 


3100 


-1.5 


3200 


-1.5 


3000 


-1.5 


3200 


-1.0 


APMPM 0559-2903 


esdM7 


3100 


-1.5 


3200 


-1.0 


3000 


-1.5 


3100 


-1.0 


LEHPM 2-59 


esdM8 


3000 


-1.5 


3100 


-1.5 


3000 


-1.5 


3100 


-1.0 








iNcar 


Intrarcd 












LP 589-7 


sd/esdM5 


320O 


-2.0 


33O0 


-1.5 


3200 


-2.0 


3200 


-1.5 


LSR 0822-1-1700 


esdM6.5 


2900 


-2.0 


3000 


-2.0 


2900 


-2.0 


3000 


-2.0 


APMPM 0559-2903 


esdM7 


3000 


-2.0 


3100 


-1.5 


3000 


-1.5 


3000 


-1.5 


LEHPM 2-59 


esdM8 


2800 


-2.0 


3000 


-2.0 


2800 


-2.0 


3000 


-2.0 



^Values in bold denote overall best fit models. 
'^Logarithm of surface gravity in cm s~^. 
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TABLE 7 
Properties of LEHPM 2-59. 



Parameter Value Rcf. 





Q4h Qg//g4 


1 




-22° 45™ 08'.'4 


1 


SpT 


esdMS 


2,3,4 


Reso'' 


18.72 


5,6 


RuSNO^ 


18.82 


7 




16.86 


5,6 


luSNO^ 


17.24 


7 


Idenis 


16.74±0.09 


8 


J 


15.52±0.05 


1 


H 


15.25±0.08 


1 


Ks 


14.76±0.11 


1 


e 


0^'746±0'.'016 yr-l 
174?7±1?2 


6 
6 


d-est 


66±15 pc 


2,4 


Vrad 


79±8 km s-i 


2 


[U, V, W]" 


[135,-180,-80] km s"! 


2 


[M/H] 


2800-3000 K 


2 


-1.5 to -2.0 


2 


References. — (1) 2MASS; (2) 



This p aper; (3) Gizi ^ 119971) : (4) 
ILepine R ich. & Shara 1 2003 d); (5) SSS 
I Hamblv et a l. 2001a b c); (6) Pokornv et alj 
J2004I) ; (7) USNO-Bl.O (Monet et al. 2003); 
(8) DENIS iEpchtein et aljil997i\ 

^2MASS coordinates, equinox J2000 and 
epoch 1998 Nov 29 (UT). 

Photographic R (IIIaF) and In (IV-N) 
magnitudes. 

Assuming a distance of 66 pc. 
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